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Abstract 

We report differential cross sections for the production of (2010) produced in 
500 GeV/c 7r~-nuclcon interactions from experiment E791 at Fcrmilab, as functions 
of Feynman-x {xp) and transverse momentum squared (pf^). We also report the D*^ 
charge asymmetry and spin-density matrix elements as functions of these variables. 
Investigation of the spin-density matrix elements shows no evidence of polarization. 
The average values of the spin alignment are [rj) = 0.01 ± 0.02 and —0.01 it 0.02 for 
leading and non-leading particles, respectively. 

Key words: PACS 13.85.Ni, 13.88.+e, 14. 40. Lb; charm hadroproduction; 
differential cross section; spin alignment; polarization 



1 Introduction 

(2010) (J^ ~ 1") production is important to understanding the produc- 
tion of charm because D*'s are expected to dominate the charm cross-section: 
their spin causes them to be favored threefold over the low-lying D mesons 
at high center of mass energies. Charm production is an interesting topic in 
its own right, being calculable in perturbative QCD[1]. There has also been 
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interest in the intrinsic charm content of hadrons[2]. Finally, accurate simu- 
lations of charm production required by future experiments at the Tevatron, 
the LHC, and lepton colliders will depend on these measurements which can 
be used to tune Pythia[3] and other Monte Carlo simulation packages. 

Since D* mesons have spin 1. we might expect the charm quark spin to be 
reflected in the meson. Spin retention is particularly strong for heavy quarks, 
where the usual stiffness of the spin vector under relativistic boosts may be 
expected to play a part [4]. In any case, there are many anomalies of spin 
and polarization. For instance, it is known that some hyperons are strongly 
polarized in some regions of production and not in others, the reasons for 
this behavior being still not completely understood [5]. It is also known that 
most of the spin of the proton is not carried by the valence quarks [6]. Copious 
production of a vector meson provides us with a unique opportunity to add 
new information on polarization, and we hope that these data will shed further 
light on hadron physics. 

In this paper, we report measurements of the production of 1^*^(2010) in 
500 GeV/c 7r~-nucleon interactions. Forward differential cross sections and 
D*^ production asymmetries have been measured as functions of Feynman-x 
{xp) and transverse momentum squared {p^)- The measurements are com- 
pared to Pythia Monte Carlo models. The differential distributions are fit to 
well-known functional forms and fit parameters are compared with previous 
measurements. Also, spin-density matrix elements were measured, again as 
functions oi xp and pf.. 



2 Experiment and Data Sample 



E791, a high statistics charm physics experiment, took data at Fermilab's 
Tagged Photon Laboratory during the 1991-2 fixed-target run. The experi- 
ment [7] used an upgraded version of the two-magnet spectrometer previously 
used in E516, E691, and E769. A 500 GeV/c tt" beam was directed at five 
target foils: a 0.52 mm thick platinum foil, followed by four 1.56 mm thick 
carbon foils, with a typical spacing of 15 mm between foil centers. Tracks and 
vertices used hits in 23 silicon microstrip and 45 wire chamber planes. The 
spectrometer included two Cerenkov counters, electromagnetic and hadronic 
calorimeters, and two muon scintillator walls. A loose trigger based on trans- 
verse calorimeter [8] energy was employed. An innovative data acquisition 
system [9] recorded 2 x 10^° events. The resulting 50 Terabyte data set was 
reconstructed using large computing farms at four sites [10]. 

The D* sample used in this analysis was collected from the decay mode D*+ — > 
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Fig. 1. D sample passing selection criteria used in this analysis. In the fit shown, 
a double Gaussian was used to model the signal shape and random event mixing 
was used to model the background. This fit yields 17800 ib 185 D*^ signal events. 

L)°7r+.Q The candidates were obtained in the decay modes K 
and D° K^n^n^n^ . Selection criteria were chosen to optimize S/\/S + B, 
where S and B are respectively the signal and background in D*"^ Q- value 
histograms. Here, the D*^ Q-value is defined as m{D*) — m{D^) — 'm{7c~^) 
where m{D*) is the reconstructed D*^ mass, m{D^) is the reconstructed 
mass, and m(7r^) is the charged pion mass [11]. Owing to the high statistics 
of the sample and the small number of criteria used for selection, the data set 
was used directly for this optimization. 

The short lifetime of the D*^ makes it impossible to reconstruct a D*^ decay 
vertex separate from the primary vertex, so in what follows, the secondary 
vertex is defined to be the point where the decays, and the primary ver- 
tex to be the point where the D*^ and are produced. The longitudinal 
separation between the primary and secondary vertices was required to be at 
least 8 (7) standard deviations^, and the secondary vertex was required to 
be outside any target foil. The component of the momentum vector of the 
candidate perpendicular to the line joining the primary and secondary vertices 
was required to be less than 0.4(0.45) GeV/c, and also the momentum vec- 
tor was required to point back to within 0.1 mm of the primary vertex in the 
plane transverse to the beam. To reduce background caused by false secondary 
vertices or tracks being mis-assigned to the secondary vertex, it was further 
required that the product over all decay tracks of the ratio of the transverse 



Charge conjugates are implied for all decay modes discussed in this paper. 
^ In some instances, different cuts were applied for the Ktt and Kutt-k modes. In 
these cases the Ktttt-k cut is shown in parentheses. 



4 



impact parameter at the secondary vertex to that at the primary vertex be 
less than one. To further reduce background, the reconstructed mass of the 
candidate was required to fall within a mass window which was determined 
by creating histograms of the masses of the D° candidates separately in the 
various regions of D*^ xp and to which a fit was performed where the 
signal was represented by a Gaussian and the background by a straight 
line. The limits of the mass window were chosen to be at 2.5<t on either side 
of the mean. Owing to the small phase space available to the D*"^ decay, the 
background in the Q- value histogram is already low enough that no additional 
cuts on the D*^ decay pion candidates were necessary, other than basic track 
quality requirements applied to all decay tracks. Applying the above selection 
criteria yields a data sample of 17800 ±185 D*"^ signal events which are shown 
in Figure 1. 



3 Differential Production Cross Sections 



For the differential cross section study, the data sample was divided into 12 
bins each of xp and p^. In the range —0.1 < xp < 0.4, the bins have a width 
of 0.05, and in the range 0.4 < xp < 0.6 the bins have a width of 0.1. In 
the range < pf. < 2 (GeV/c)^ bins have a width of 0.5 (GeV/c)^, and 
in the range 2 < p"^ < 10 (GeV/c)^ bins have a width of 1 (GeV/c)^. In 
each bin of xp or D*"^ Q- value histograms were created. Random event 
mixing was used to model the background shape, and the background level 
was determined by fitting the event mixed histogram to the corresponding 
data histogram using a binned maximum likelihood fit in the background 
region (10 < Q < 20 MeV/c^). In the signal region (0 < Q < 10 MeV/c^) we 
subtract from the data histogram the scaled number of events from the event- 
mixed histogram to obtain the number of signal events. In this way, it isn't 
necessary to model the details of the signal shape to determine the number of 
signal events in the histogram. 

To model the acceptance, 14 million Monte Carlo events were generated which, 
when reconstructed and subjected to the same selection criteria as our data, 
yielded approximately ten times as many D*"^ events as in the data sample. 
The Monte Carlo samples were fit in the same way as the data to obtain a 
bin-by-bin efficiency estimate. 

The acceptance-corrected yields were used to plot the differential cross sections 
in Figure 2. For the xp distributions, we used data in the range < p^ < 
10 GeV^, and for the pf. distributions, —0.1 < xp < 0.6. The xp distribution 
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Fig. 2. Acceptance-corrected differential cross sections of D*^ production as func- 
tions of (left) and (right). Error bars include both statistical and systematic 
errors. Fits shown are of functional forms described in Equations 1 and 3. Also, 
default and -D^-tuned Pythia models are plotted. 

was fit in the range 0.1-0.6 with the traditional function: 



da 



(1) 



The distribution was fit with the traditional function 

da 



dp 



,2 



exp(-5p^) (2) 



in the range < pf. < 4 (GeV/c)^ and with the function 

da ( 1 
dp^ \bmc^ + 



da ( 1 \^ . . 

2 ^ I I, ■}. I „2 ) V"^'' 



in the full range < pf. < 10 (GeV/c)^. This last functional form was sug- 
gested by Frixione et al.[12]. 

In Figure 2, the E791 results are also compared to the default Pythia model [3] 

and a "tuned" Pythia model [13] in which parameters are adjusted to more 
closely model the D"^ production characteristics reported in an earlier E791 
paper [14]. The results of the various fits are listed in Table 1, along with 
results from previous experiments [15-19]. 



6 



Tabic 1 

Parameters from fits to tt^N differential D*^ cross sections vs. xp and pj^, including 
comparisons to other experiments ([15-19]). Column labels refer to fit parameters 
of the same name defined by Equations 1-3. The labels "1" and "nl" refer to leading 
and non-leading particles, respectively. Asterisks (*) indicate extended range of 
used in fit of Equation 3. Errors shown on E791 results include both statistical and 
systematic components. 



Expt. 


P 

beam 

GeV/c 


Xp fit 
range 


n 


fit range 
(Gev/c)2 


B 


hmc? * 
(GeV/c)2 


/? * 


E791 


500 


0.1 to 0.6 


3.8 ±0.2 


to 4 


0.75 ± 0.02 


5.1 ±0.8 


5.0 ±0.6 


1. (D*-) 






3.6 ±0.2 


or 


0.73 ± 0.02 


6.1 ± 1.5 


5.6 ± 1.0 


nl. {D*+) 






4.1 ±0.3 


to 10* 


0.77 ±0.03 


4.0 ± 1.1 


4.4 ±0.8 


WA92[15] 


350 


to 0.6 


4.3 ±0.4 


to 4 


0.84 ± 0.05 


3.6 ±0.8 


4.7 ±0.7 


1. 






4.9 ±0.5 


or 








nl. 






3.9 ±0.4 


to 14* 








E769[16] 


250 


0.1 to 0.6 


3.5 ±0.3 


to 4 


0.70 ±0.07 






1. 






2.9 ±0.4 




0.58 ± 0.09 






nl. 






4.1 ± 0.5 




0.79 ± 0.09 






NA32[17] 
1. 

nl. 


230 


to 0.8 


o 14+0.40 

9 /^o+0.53 
^•0^-0.49 

Q OQ + 0.66 

'J-O'J-0.62 


to 10 


0.79 ± 0.07 

Q 71+0.09 
'1-0.08 

0.90 ±0.11 






NA32[18] 
1. 

nl. 


200 


to 0.7 


2-8lo!9 

4 7+1-9 
^•'-1.6 

1 7+1-^ 
-■-•'-1.0 


to 5 


u.y_o.2 






NA27[19] 


360 


to 0.5 




to 3 


0.9 ±0.4 







4 Chcirge Asymmetry in Production 



The charge production asymmetry is defined by the parameter: 

_ (7(D*-) - a{D*+) _ N{D*-) - N{D*+) 

" a{D*-) + a{D*+) ~ N{D*-) + N{D*+) ^ 



where a{D*^) and a{D*^) denote the production cross-sections, of D*~ and 
respectively, and N{D*~) and N{D*~^) denote the respective acceptance- 
corrected numbers of such particles observed. 
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Fig. 3. Charge production asymmetry as functions oixp and p\-. (Error bars include 
statistical and systematic errors added in quadrature. Horizontal error bars indicate 
RMS values of Xi? or distributions within the bins.) Comparisons to default and 
D^-tuned Pythia D*"^ asymmetry predictions also shown. 

Using the same fitting techniques and data sample as was used in the differ- 
ential cross section study, the production asymmetry as functions of xp and 
are shown in Figure 3. 



5 Spin-Density Matrix Elements 



The polarization state of a spin 1 particle can be described by a complex 
3x3 spin-density matrix. A complex 3x3 matrix has 18 real components, 
but hermiticity and the fact that Trp = 1 reduces this number to 8. Also, the 
D*^'s are produced in a strong interaction where parity conservation further 
reduces this number to 4. The spin-density matrix when expressed in a helicity 
basis then takes this form [20] : 



Pii Pio Pi-i 
Pw 1 - 2pii -p\q 
yPi-i —Pio Pll J 



with and of course pn real. 

Prom this, one can derive the angular distribution [20] , 
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^._Lli + (i_3oos^9)(3pn-l 

sin 29 cos0)9?pio 
— (3sin^ 6'cos20)pi_i]. 



(5) 



The angles 6 and (6 in the above equation are polar and azimutlial angles 
respectively in a right-handed spherical polar coordinate system in the rest 
frame of the D*"^, whose polar axis {9 — 0) points along the D* boost vector 
(a "helicity" basis, with the z-axis chosen to lie in the production plane) and 
the (f) — 0,6 = 7r/2 axis is chosen to point in the direction d x where d is a 
unit vector in the direction of the D* momentum vector and la is a vector in 
the direction of the beam particle momentum vector (the x-axis is normal to 
the production plane). 

Note that this distribution depends on only three of the four independent 
components of the spin-density matrix: pn, 3?pio, and Thus Q'pio cannot 
be determined from this angular distribution study. 

The degree of polarization is commonly expressed by the spin alignment pa- 
rameter: 

77 = 1 - 3pii 



To measure the spin-density matrix elements, the data sample was first divided 
into five bins each of xp and p^. In each of these bins, an 8 x 8 two-dimensional 
histogram of da/dQ vs. (cos^, 0) was created. The D*^ yield in each bin was 
measured as follows. Events were counted both above and below a Q-value of 
10 MeV/c^ for both the data histogram and the event-mixed data histogram 
(our background model). A comparison of the number of events in the data 
and mixed data sample in the > 10 MeV/c^ region yields a scale factor, which 
is used in subtracting the number of background events from the signal region 
< 10 MeV/c^. The same procedure was applied to a Monte Carlo sample to 
obtain a binned acceptance function, and the data sample was corrected for 
acceptance in this way. 

The resulting two dimensional histogram was fitted with the functional form 
shown in Equation 5 to obtain values for the spin-density matrix elements. 
These results are plotted separately for D*~^ and D*~ in Figure 4. Because 
the sample had to be split into sub-regions of cos 9 and (p in order to perform 
the fits, fewer bins in xp and were used than in previous sections. There is 
no measurable polarization at our level of sensitivity in any sub-region of the 
xp or range examined. The average value of the spin alignment parameter 
is found to be {r]) = 0.01 ± 0.02 and -0.01 ± 0.02 for leading (D*-) and 
non- leading (/?*+) particles respectively. 
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Fig. 4. Spin-density matrix elements as functions of a;^ and p^. (Error bars in- 
clude systematic as well as statistical errors added in quadrature.) The solid line 
corresponds to unpolarized D*^s. 



6 Systematic Errors 



The systematic errors reported were determined by comparing results, either 
from various subsamples, or by analyzing the full sample in multiple ways. 
Differences in measurements of the number of observed D*^ signal events 
were then determined as fractions of statistical errors bin by bin. Root mean 
squares of these fractions were then used as estimates of the systematic error 
as a fraction of statistical error. Several possible sources of systematic error 
were considered, and those found to be significant are tabulated in Table 2. 
The cut on the mass varies with xp and p'^ and depends on our modeling of 
the signal. By comparing results with different mass cuts, a systematic 
error of 24% (of the statistical error) was assigned for this effect. Because the 
default Pythia and xp distributions don't match those in data exactly, it 
was necessary to weight the Monte Carlo events to correct the distributions 
when modeling acceptance as a function of the other variable. By varying the 
"weighting curve" parameters by la from their central values, it was deter- 
mined that a relative systematic error of 13% should be assigned. Two decay 
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Table 2 

Summary of systematic errors reported as fractions of statistical errors. 



Source of Uncertainty Syst. / stat. err. 

Pt/xf weighting 13% 

£)° mass cut 24% 

Tracking efficiency 25% 

D* fitting 39% 

Drift chamber efficiency modeling 58% 

Total 79% 



modes of the (Ktt and KinrTr) are used in this analysis. Since one mode in- 
volves two more decay particles than the other, one can compare results from 
the two modes to obtain an estimate of the error resulting from modeling of 
the tracking efficiency. A relative error of 25% was assigned for this effect, 
although some uncertainty remains as the two samples were statistically inde- 
pendent, and it is not possible to unambiguously untangle the systematic and 
statistical errors. Event mixing was used to model the backgrounds in the D*^ 
histograms. If the event-mixed histograms do not perfectly model the actual 
background, this would result in a systematic error. By varying the range of 
Q- values in which the background is fit, an estimate of the magnitude of the 
systematic error is obtained. A relative error of 39% is assigned to this effect. 
Over the period of E791's run, the 500 GeV/c vr" beam caused a significant 
amount of ionization in regions of the drift chambers through which it passed. 
This caused deposits to gradually form on the drift chamber wires resulting in 
a loss of tracking efficiency in regions through which the beam passed. These 
drift chamber "holes" changed in size and shape throughout the run period. 
Although this is modeled in the Monte Carlo in a time dependent fashion, 
this modeling is not perfect. By comparing results obtained from subsamples 
of various run periods, a 58% relative systematic error was assigned to this 
effect. 

Each source of error was measured as a fraction of the statistical error, and a 
quadratic sum of the various sources is shown at the end of the table to be 79% 
of the statistical error. Combining the statistical with the overall systematic 
error by adding them in quadrature results in a combined error which is 1.3 
times the statistical error alone. These combined errors are used to create the 
error bars in the plots in this paper. 
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7 Conclusions 



We have measured the differential cross sections of D*^ production as func- 
tions oixp and p^. While the usual functional forms used by other experiments 
to parameterize their data (Equations 1 and 2) do not describe our data well, 
Equation 3 suggested by Frixionc ct al. provides an excellent fit over our full 
range of (x^ = 8.4 for 9 degrees of freedom). The Pythia model tuned for 
the E791 analysis of production does describe the cross section much 
better than the default Pythia model. 

We observe an overall charge asymmetry which favors the leading particle as 
predicted by the Pythia model, although our overall observed asymmetry is 
significantly less than predicted. We do observe an increase in asymmetry with 
xp as predicted by Pythia and consistent with that for production [14]. 
The asymmetry is approximately constant as a function of while it dips 
to zero near = 0, data points in this region are also consistent with the 
average value. 

We observe no D*"^ polarization at our level of sensitivity, either on average 
or as functions oi Xp or p^. 
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